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9.1 Introduction 



In a thermoelectric material, heat can be transported or used to generate electricity based on the Peltier and 
Seebeck effects. The efficiency of a thermoelectric device is traditionally described in terms of the extensive' 1 
or system parameters such as hot- and cold-side temperature, length and area of thermoelectric element, 
and applied voltage or load resistance. In only the most simplified cases (e.g., temperature-independent 
thermoelectric properties) can this efficiency be computed analytically, where it can be shown that the 
thermoelectric figure-of-merit z is the intensive material property of prime importance (see Chapter 1, 2, 
and ll). 1 ' 2 At the optimal electric current, z alone determines the efficiency. 

Although very instructive, such simplifications ignore the effect of thermoelectric compatibility in real 
thermoelectric devices. The thermoelectric compatibility factor is the reduced electric current, which is 



'Extensive refers to properties that depend on material sample size, 
intensive refers to properties independent of material sample size. 
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necessary to achieve the highest efficiency determined by z. Because the compatibility factor changes with 
temperature while the electric current is constrained, the efficiency of a real device will be less than that 
calculated from z. The effect of thermoelectric compatibility is most important for segmented 
thermoelectric generators, but it also affects the exact calculation of performance for all thermoelectric 
devices. 

To calculate the exact performance of a thermoelectric generator analytically, it is simplest to use a 
reduced variables approach that will separate the intensive properties and variables (such as temperature 
gradient, Seebeck coefficient, current density, heat flux density) from the extensive ones (e.g., voltage, 
temperature difference, power output, area, length, resistance, load resistance)/' This approach allows a 
definition of a local, intensive efficiency in addition to the traditional system efficiency 5-8 as well as the 
derivation of the compatibility factor. 4 

9.1.1 Definitions 

Consider the one dimensional (V = — ), steady-state, thermoelectric power generation problem, where 

ax 

(for now) only a single (n- or p-type) leg is considered. The thermoelectric material properties all vary 
with (absolute) temperature T: the Seebeck coefficient a, the thermal conductivity k, and electric 
resistivity p. In the following discussion, we will assume isotropic materials properties. For anisotropic 
materials, the tensor representations of the materials properties and equations are required. Positive 
electric current density/ > and heat flux (heat current density) Q > (with units of Watt/cm ) flows 
from Tj, to T c (Figure 9.1). Positive electric field E and temperature gradient VTare in the opposite 
direction of/ and Q. The subscripts h and c (lower case) denote the value at a particular (hot- or cold- 
side) temperature (K h = K(T h ), VT h = VT(T h ). 

The electric current density is for a simple generator, given by 



J 



1 
A 



(9.1) 



where / is the electric current (Amps) and A is the cross-sectional area of the thermoelectric 
element. 

The electric field is given by a combination of the reversible Seebeck effect and the 
irreversible effect of Ohm's law. Using the sign convention described above, the electric field 
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FIGURE 9.1 Diagram of a single-element thermoelectric generator. The direction of positive variables is shown 
relative to the hot- and cold-side. For positive Seebeck coefficient (a), all of the variables are positive for a generator 
operating efficiently. For negative Seebeck coefficient (a < 0), the electric current, field, and potential (/, E, V) will be 
negative, or opposite to the direction shown. 
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from a purely resistive element using Ohm's law is E — —pj. The electric field produced by the 
Seebeck effect is E — aVT. Combining the Seebeck and Ohm effects, gives the electric field at any 
position: 

E = aVT - pj (9.2) 

Similarly, heat is transported reversibly by the Peltier effect, Q = aTJ (where aT is the Peltier 
coefficient, related to a by the Onsager reciprocal relations ) and irreversibly by Fourier's law 
Q = kVX (using sign convention of Figure 9.1): 

Q = aT] + kVT (9.3) 

The Peltier effect is often considered a surface effect between two materials but the heat transported 
is a property of a single material. 

In both cases (Equation 9.2 and Equation 9.3), the irreversible and reversible effects are treated 
independently, and can simply be summed. This treatment is related to Kelvin's assumption. The 
irreversible heat flow is further constrained by the steady-state heat equation 

V(kVT) = -T-^JVT - pj 2 (9.4) 

dr 

da 
where T- — is the Thomson coefficient. From Equation 9.4, we see that the heat produced (or 

consumed) by the Thomson effect is transported away from (or toward) the location by the irreversible 
part of Equation 9.3 as a part of kVT. 

The electric power density P (power produced per volume) is the product of the electric field E and 
current density /: 

P = EJ (9.5) 

Using the sign convention in Figure 9.1, a purely resistive element (a = 0) would require a negative 
electric field E = —pj to make a positive current (+/) so that the power density P = — pj 2 is negative 
(electric energy consumed). 

The heat equation (Equation 9.4) which includes the Thomson effect can be derived by invoking the 
conservation of energy. The divergence of the heat flow should be compensated by the sources (or sinks) 
from the electric power generation (work): 

VQ = P (9.6) 

Substituting Equation 9.3, Equation 9.2, and Equation 9.5 into Equation 9.6 and evaluating V/ and 
V(aT) will derive the steady-state heat Equation 9.4. For the one-dimensional problem 

V7 = (9.7) 

because there is no buildup of electric current in the steady-state. To evaluate V(aT), it is first 

noted that, in the power generation problem, there is a temperature drop at every point (VT > 0, 

for thermoelectric material and electric and thermal contacts). Thus, the temperature profile T(x) is 

invertible to a well defined x(T). In this way, a(T) = a(x(T), T) is defined and then the gradient of 

da 
a is simply Va — - — VT (note the use of total not partial derivatives). This gives 
dT 

V(aT) = T-^-VT+aVT (9.8) 

dr 
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9.2 Reduced Variables 
9.2.1 Relative Current Density 

For calculating and comparing generator properties it is most convenient to use a reduced variable in 
place of the actual current density /. This is because the most efficient current density scales with the 
length of the thermoelectric element in the same way as VT does. By dividing the two, a simplified 
expression for reduced efficiency is derived. 

Here, we use is the relative (reduced) current density u, which is the ratio of the electric current density 
(/) to the heat flux by conduction (kVT) given by 

(9.9) 



kVT 



(again, VT > for power generation). The electric power density (from Equation 9.2 and Equation 9.5) 
in terms of u is given by 

P = k(VT) 2 X u(a - upn) (9.10) 

The heat flux density (from Equation 9.3) in reduced variables is 

Q = kVT(cmT+\) (9.11) 

The variation of u is governed by the heat Equation 9.4. For the one- dimensional problem with cross- 
da 

sectional area A a constant, Equation 9.4 combined with Equation 9.7 gives V(1/m) = — T— — VT — pj. 

dT 

For this one-dimensional problem we can substitute V(l/u) = — r-r~' ~T~ = "1 — VT an d J = ukVT 

u ax ax dT 

to give 

du 7m da , 

= u T h u ok (9.12) 

dT dT ^ ' 

Again, the material properties a, k, p can be functions of temperature. Notice that here we have removed all 
reference to the spatial coordinate x so that u can be considered a function of temperature Tonly. 

9.2.1.1 Other Representations of the Reduced Current 

/ . 3 

The relative current density u = is the most instructive form of the reduced current. Sherman uses 

K\T 

y ~ 1/u, which simplifies the thermoelectric potential (Equation 9.21) but becomes ill defined for an 
open circuit generator when / = u = 0. For power generation, small u is of most interest, specifically 

< u ^ zla (Figure 9.2). Other multiplicative factors of — — are less instructive because u has the least 

variation in a thermoelectric generator and therefore allows the comparison of compatibility factors 

( du \ 
(Equation 9.18). For example, the relative change I — - — ) of u is vanishingly small for small u 

\udT } 

du I da \ . i 8 

(7 = 0, u ~ 0): from Equation 9.12, — - — = u\ \- upK I vanishes as u approaches zero. Other forms 

udT \ dT ) 

a it/ 

of the reduced current such as u- or i= — = aTu have nonvanishing relative derivatives, 

di da I J da i \ • , • 

e.g., = 1 1- il 1 r- 1 does not vanish as i approaches zero. 

z'dT adT T \ adT zT l ) 
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FIGURE 9.2 Variation of reduced efficiency (Equation 9.15) with relative current density, u. The maximum 
efficiency is achieved at the compatibility factor, u = s. For the plot, zT = 1, aT = 0.1 V similar to the values for 
(Bi,Sb) 2 Te 3 . 

9.2.2 Reduced Efficiency 

Efficiency 17 is defined as the power produced divided by the power supplied to the system. The 
(infinitesimal) efficiency along the infinitesimal distance Ax is the power produced (per cross-sectional 
area) Pdx divided by the heat flux through, Q, or 

Pdx 

~Q 

dr 



17 



In Equation 9.13, dx is in the direction of the temperature gradient so dx 
Equation 9.10 and Equation 9.11 into Equation 9.13 gives 

dT u(a — upK) 
1= — f— 



VT 



(9.13) 
Substituting this with 



T 



(9.14) 



T 



AT 



The first term is recognizable as the infinitesimal Carnot efficiency, rj c = . The reduced efficiency, 

17,., defined by rj = ri c rj r , is not an infinitesimal quantity: 



u(a — upK) 
~1~ 



(9.15) 



T 



The reduced efficiency can be succinctly written in terms of u[ — I when z^0, where z is the 
thermoelectric figure-of-merit: 



If the current is nonzero, u =£ 0, (/ 5^ 0), and a ¥= 0, the reduced efficiency is simply written as 

a 



(9.16) 



1 - u- 



i?r 



1 + 



1 



(9.17) 



uaT 



Similar to u, the reduced efficiency r) T (u(T), T) is simply a one-dimensional function of temperature once 
an initial u is applied. This allows the exact calculation of efficiency using a simple spreadsheet calculation 
(described in Section 9.4.2). 
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9.2.3 Efficiency Dependence on Current 

Whether in power generation or Peltier cooling mode, the reversible, useful thermoelectric effects 
compete with the irreversible Joule heating. Because the linear effects are directly proportional to the 
electric current while the irreversible Joule heating is proportional to the square of the current, there is 
necessarily an optimum operating current to achieve the optimum efficiency. The variation of reduced 
efficiency with u current (Figure 9.2, Equation 9.15) is analogues to the variation of the power output to 
the electric current: at zero u current, there is voltage produced but neither power nor efficiency. As u 

increases, the efficiency increases to a maximum value and then decreases through zero. Past this zero- 

z 
efficiency crossing where u = — , the Ohmic voltage drop is greater than the Seebeck voltage produced, 

a 

and thus the power output and efficiency are negative. 

The value of u which gives the largest reduced efficiency (Equation 9.15) is thermoelectric (power 
generation) compatibility factor s: 

•Jl+zT - 1 

(9.18) 



aT 



For small zT, this can be approximated by 



^- (9.19) 

2a 



This largest reduced efficiency r) t (u — s) is given by 



Vl + zT - 1 

maX1 ?r = 7= S— (9.20) 

VI + zT + 1 

Thus, in the most general case (a, k, p, temperature dependent), the thermoelectric figure-of-merit z is 
the material property that determines the maximum local efficiency. This can be used to derive the 
definition of z as the local thermoelectric figure-of-merit. 

From Equation 9.18 it is clear that the compatibility factor s is, like z, a temperature-dependent 
materials property derived from the temperature-dependent materials properties a, k, p. Thus s cannot 
be changed with device geometry or the alteration of electric or thermal currents. 

If u t^ s then the efficiency is less than the maximum efficiency of Equation 9.20. Since 

u = , there is some control over u from the applied current density / (traditionally regulated 

kvT 

by a load resistance). However, once u is selected at one point, it cannot be adjusted in a 
thermoelectric element to follow the temperature variation of s (Figure 9.3), because the variation of 
u is fixed by the heat equation (Equation 9.12). 

Conveniently, the variation of u within a thermoelectric leg is typically small. Since all segments in a 
thermoelectric element are electrically and thermally in series, the same current I — A] and similar 
conduction heat AkVT flow through each segment. When the electric current is near zero (J ~ 0) the 
heat flow is very uniform (V(kVT) ~ 0) so u is nearly constant. For 7 t 4 0, the conduction heat is only 
slightly modified by the change in temperature gradient due to the Thomson and Joule sources of heat 
(Equation 9.4). Thermoelectric generators operating at peak efficiency typically have u that varies less 
than 20% within all thermoelectric materials in the entire element (Figure 9.3). 

To a reasonable approximation, u, once established, remains constant throughout the thermoelectric 
element. 1 Thermoelectric coolers, on the other hand, are typically driven with much higher w, up to 
u = oo (VT = 0) at the end being actively cooled. 

The actual reduced efficiency of a material depends not only on the maximum reduced efficiency 
(Equation 9.20) determined by z, but also on how close u is to s (Figure 9.2). The actual reduced efficiency 
(Equation 9.15) is always less than the maximum reduced efficiency (Equation 9.20), because u, as 
determined by the heat equation (Equation 9.12), varies differently from the material property s, so they 
cannot be equal at more than a few isolated points. The difference between the maximum and actual 
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FIGURE 9.3 Variation of relative current density u with temperature for a typical thermoelectric generator. The 
total variation of u within a material and the change at the segment interfaces is less than 20%. The u shown is that 
which gives the highest overall efficiency. For (Bi,Sb) 2 Te3 and Zn 4 Sb 3) u is less than the compatibility factor s, while for 
the CeFe 4 Sb 12 segment u is greater than s. 



reduced efficiency is largest for large differences between u and s. This can be seen graphically in Figure 9.3 
and Figure 9.4. 



9.2.4 Thermoelectric Potential 

For many expressions and calculations, it is convenient to use the thermoelectric potential <P (with units 
of Volts), 4 which is related to the electrochemical potential (jl (pb = q<P with q the charge/particle). The 
thermoelectric potential is given by 



<P = aT + \lu = aT 



kVT 



The heat flux Q containing both the Peltier and Fourier heat (Equation 9.3) is simply 

Q = J<P 



(9.21) 



(9.22) 
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FIGURE 9.4 Local, reduced efficiency (using optimized u from Figure 9.3) compared to the maximum reduced 
efficiency (if u = s for all temperatures, Equation 9.20). The difference is most substantial in the regions where u is 
most distant from s (Figure 9.3). 
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Multiplying by A gives the heat flow U — AQ in terms of current: 

U = I<P (9.23) 

The electric field E containing both the Seebeck and Ohm effects (Equation 9.2) is given by 

E = V<J> (9.24) 

where Equation 9.4 and Equation 9.8 are used. Thus, the electric power density (Equation 9.5) 
produced is 

P = 7V0 (9.25) 



9.3 Generator Efficiency and Performance 
9.3.1 Efficiency of a Thermoelectric Device 

Once the intensive (local) reduced efficiency is known at every point, the total system efficiency can be 
calculated. Because efficiency is strictly a fractional quantity less than one, care must be taken to use the 
appropriate summation metric when calculating combined efficiencies. 

For a general energy conversion process, let U be the heat energy into the system and W be the work 

removed. By conservation of energy, the heat energy out of the system is [/,- — W, 17, = — - . The 

u i 
efficiency of such a system is 

W 
V = — (9.26) 

9.3.1.1 Parallel Efficiency 

For two thermoelectric elements thermally in parallel (such as a thermocouple of n- and p-type 
elements), the combined efficiency is a (weighted) average of the efficiency of both generators, weighted 
by the heat flowing through each generator: 

_ W r + W 2 _ Villi + V2U2 , Q __ 



9.3.2 Series Efficiency 

For two processes in series the heat out of system one is the heat into system two: [7 2 = U\ — Wy. The 

efficiency of the combined system is the work from both systems divided by the energy supplied to the 

W l + W 2 
first system in the series 171+2 series = • Combining these relations gives the summation rule for 

efficiencies in series: 

1 " T?l+2,series = (1 ~ 1?l)(l ~ Ife) (9.28) 

For many processes in series this becomes 

1 " Series = f[ l " ^ ^ 9 ' 29 ) 

i 

By taking the logarithm of both sides, the product series can be changed to a summation, which 
in the limit of infinitesimally small steps becomes an integral. For infinitesimally small rj h we 
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have ln(l — 17,) = — 77,. ' 



ln(l - T7 series ) = I ^ ln(l - 17;) = - I 77; (9.30) 

Note that since i7 series 7^ J, ,17,, it is best to refrain from identifying infinitesinially small 17, with drj; 
instead, it is simply a local efficiency: 



' - 1?local 



^series = 1 ~ exp - 17 local (9.31) 

Using the expression for the local efficiency in terms of the reduced efficiency (Equation 9.14), the overall 
efficiency of a finite segment can be derived using temperature T to define the integration path: 



H>] 



rj=l-exp|- -^-dTl (9.32) 



9.3.3 Single Thermoelectric Element 

A thermoelectric generator consists of a n-type element and p-type element connected thermally in 
parallel but electrically in series. Often, it is simpler to compare the efficiency and performance of 
individual elements rather than n-p couples, for example, when one is selecting materials. 

9.3.3.1 Thermoelectric Element Efficiency 

The thermoelectric potential becomes a convenient integration variable for evaluating Equation 9.31. 
Integrating the local efficiency (Equation 9.13) by substituting Equation 9.22 and Equation 9.25 becomes 
simply a function of the thermoelectric potential at the two ends of the process: 

J Tc Q dx = J Tc ~W dx = J Tc -* = J Tc dl "^ = ^^ = In(^) (9.33) 

Here, the thermoelectric potential is a parametric function of T, because u is also a function of T (e.g., 
<$> c = <P{n{T c ),T c )). This provides a simple expression for the efficiency (Equation 9.31) of a single 
thermoelectric element: 



v = i - -r- = -r- (9-34) 



In terms of u and T this is 







a c T c 


1 
+ — 


«h^ h 


1 

H 

"h 



V = 1 - ^ (9.35) 



The maximum efficiency of a thermoelectric element is computed by finding the initial u that maximizes 
Equation 9.35 for a given temperature range (T h and T c ). Once an initial u (for example, M h ) is 
established, all the other u(T) (for example, u c ) are defined by Equation 9.12. 

9.3.3.2 Thermoelectric Element Performance 

d<S 

With E - V<5 = — (Equation 9.24), the output voltage V — J Edx = J d$ = <P h - <P C - A<P is 

V = A® (9.36) 
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Because u = —= — , any current density / can lead to any u provided the appropriate temperature 

gradient is supplied. The absolute magnitudes of/ and kVX, for a given u, will depend on the length 

of the thermoelectric element I — x(T b ) — x(T c ) (Figure 9.1). Integrating u= —= — , gives 

kvT 

J" ku dT = \ ] dx. Assuming a uniform element with constant cross-sectional area (J(x) = constant) 
gives: 

ku dT = Jl (9.37) 



i; 



Thus the heat flux Q (at any point) (Equation 9.22) can be rewritten 



1 ^T 1 



Q = - x ku dT (9.38) 

I J T c 

Similarly, the power output (work) per unit area I — = J P dx I is the product of Q h (Equation 9.38) 
and efficiency (Equation 9.34) can be written 

W 1 f r '> 

— = X A0 km dT (9.39) 

A Z Jr c 

The design condition of maximum thermoelectric efficiency sets the value of u(T), and therefore also <P 
and A$. Thus the heat flux and power output at maximum efficiency conditions are directly 
proportional to III, and can be computed from Equation 9.38 and Equation 9.39 once / is selected. 
Typically, / is a design variable set to work with a desired heat flux (see Section 9.6.2). 
The temperature variation along the length l(T) can be calculated from Equation 9.37: 

i r r 

l(T) = - uk dT (9.40) 

J J T c 

9.3.4 Thermoelectric Couple 

The performance and output of the entire generator scales in a simple way with that of the thermoelectric 
couple. The performance of the couple, however is not simply the average or sum of the corresponding 
n- and p-elements. 

9.3.4.1 Thermoelectric Couple Efficiency 

For an n- and p-element in parallel (Equation 9.27) the efficiency of each segment (Equation 9.34) can be 
combined using the heat flow Equation 9.23. Since is the electric current J flows in the opposite direction 
in the n-element compared to the p-element (Figure 9.6), we have I n — — L: 

(9.41) 
In terms of u and T this is 



A0 p 


-A0 n 




Vn+p q 


-#n 




ap.c^c H 

"p,c 


"n,c^ c 


1 

u n,c 


1 

a p ,h7h H 

M p,h 


- a n ,hTh - 


1 



ij = 1 - P 7- (9.42) 



For an n-type material operating efficiently, a n , u n , <P n , and A0 n are all negative. Finding the maximum 
efficiency of Equation 9.42 requires optimizing the initial conditions for both n- and p-type elements 
(for example, u pil and M nh ). 
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9.3.4.2 Thermoelectric Couple Performance 

The output voltage is the sum of the voltage from the two elements (Equation 9.36): 

V = A0 p - A0 n (9.43) 

For a typical n, p -thermocouple, the electric current through both elements is the same, but the relative 
contribution of heat flux is regulated by having different cross-sectional areas A„ and A n . The maximum 
efficiency condition, which defines uJT) and u n (T), requires a specific ratio of A p /A n . 

Assuming the n-element and p-element have the same total length, the ratio of the cross-sectional 
areas can be calculated by equating the currents (±1 = JA): 

I = J V A V = -J n A n (9.44) 

Combining this equation with Equation 9.37 gives 

f Th 
Ap - h ~ )T < - (9.45) 



J f h 



P dr 



Incidentally, the most efficient area ratio is found from the most efficient uJT) and w n (T). 

Equation 9.44 and Equation 9.45 with Equation 9.37 can be solved for /, in terms of the total area 

A to ,ai = A p +A n (9.46) 

to length ratio (regardless of whether the most efficient u are used): 



MpKpdTX u n K n dT 
u p K p dT - u n K n dT 



= ^X rT h ' 



Notice that once a u is selected the current is directly proportional to the A total /l ratio. 

The load resistance i?Load is traditionally used to adjust the current. For given m, the load resistance can 
be calculated using the current from Equation 9.47, combined with Ohm's law and Equation 9.43 to give 



A0„ - A0 n 

l Load 



^Load = ~ P T '" (9.48) 



The current can be used to find the total heat flux and power produced by the thermocouple. The 
combined, total heat flowing L/ tota i = Q to talAotal is the sum of the heat flowing through each (n- and 
p-type) element (U p + U n = I& p - /$„), giving: 

[/total = J(*p - #„) (9-49) 

Using Equation 9.26, Equation 9.41, and Equation 9.49 (or W — IV and Equation 9.43) the total electric 
power (work/time) is simply 

W = 7(A0 p - A0 n ) (9.50) 

The temperature variation along the length 1{T), which will have different variation for n- and 
p-elements, can be calculated from Equation 9.40. 
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9.4 Computation of Generator Performance 
9.4.1 Analytic Example Using Constant Coefficients 

When a, k, p, and therefore z, are constant with respect to temperature the performance of a generator 
operating at maximum efficiency can be calculated analytically. The solution to Equation 9.12 for u(T) in 
this case is 



1 



1 



2(T - T c ) K p 



(9.51) 



where u(T c ) = u c is used as the (not necessarily most efficient) initial condition. 

With T — TJ, in Equation 9.51, the efficiency (Equation 9.35) can be maximized with respect to the 
initial u condition (e.g., u^) to find the most efficient u(T): 

2 



1 



1 



s(T) 2 



(T - T)2Kp ■ 



(AT -V 

{-K P sm) 



(9.52) 



where T = 



Tu + T c 



Note that when (a, p, k) are constant with respect to temperature, z is also a 



constant but s is not (Figure 9.5). 

For w h and u c this gives — = — — 
, u h s(T) 

given by 



AT - 1 1 AT _ , , . „. 

Kps(T), — = — = — I Kps(T). The resulting efficiency is 

2 u c s(T) 2 



1 



AT Vl+zT- 1 



(9.53) 



zT+r c /T h 

This equation is normally derived starting with the extensive expression for efficiency. 

For comparison, if the maximum reduced efficiency could be used (Equation 9.20), known as an 
infinitely staged or cascaded system, where u — s, the efficiency is given by 



r,= l 



(VT 
(VT 



z t c + iy 



zT h 



If 6XP l 



2(yi + zT h - yrr^) • 

(VI + zT b + l)(Vl + zT c + 1) , 



(9.54) 



Equation 9.53 demonstrates the importance of the figure-of-merit when the compatibility factor is 
nearly temperature independent. For real materials where a, k, p, vary with temperature, it is 
commonly desired to find an averaged Z (upper case) ' ~ to calculate the efficiency in place of z 



> 3 

Si 

■D 

§ 2 
to 




/ 

z/for highest efficiency 
z= constant 

a =200 /iV7K 
p = 2 mil cm 
X-=10mW/cmK 



s= compatibility factor 



300 400 500 600 700 800 900 1000 
Temperature (K) 

FIGURE 9.5 Variation of relative current density, u compared to the compatibility factor 5 with temperature for a 
thermoelectric generator with a, k, p, and therefore z, constant. 
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(lower case) defined by Equation 9.16. It is the effect of the compatibility factor that explains why 
averaging works well in some cases (when s does not vary significantly) and not in others (when s 
varies by more than a factor of two). 

9.4.2 Calculation of Exact Solution 

The methods typically used for the computation of efficiency are complex, requiring finite element 
methods that include both volume and surface terms using averaged material parameters. ,5 ' 14 ~ 17 By 
using the reduced variables, however, computation and optimization of segmented thermoelectric 
generator performance (including even contact resistances) can be performed with a simple spreadsheet- 
type calculation. Because u, once established, can be represented by a one-dimensional function of T, the 
generator efficiency can be calculated with a one-dimensional function. 

The thermoelectric properties (a, p, k) are tabulated as a function of temperature. For now, we need 
only the temperature drop desired across a material, not the physical length. The length will be calculated 
later and the temperature drop allocated can be adjusted to coincide with a desired length if necessary. In 
Table 9.1 and Table 9.2 the properties are tabulated at a minimum of every 25 K. 

For computation, the differential Equation 9.12 can be approximated by combining the zero Thomson 
effect (daldT = 0) solution with the zero resistance (pK = 0) solution 1 ' : 



1 1 

U„ M„_i 



y]l - 2u 2 n _ lP J<AT - TAa (9.55) 



where Aa = a(T n ) — a(T„_j) and p~R denotes the average of pK between T„ and T n _ l . 

At the interface between two materials, where a may be discontinuous, the properties of both materials 
are tabulated at the same interface temperature (AT = 0, Aa ^ 0). This allows the discontinuous change 
in u to be correctly calculated from Equation 9.12 (Aw = u 2 TAa). 

Using Equation 9.55, u{T) can be calculated given an initial condition (« h is used in example of 
Table 9.1). The maximum single-element efficiency is found by varying these initial u conditions and 
calculating the efficiency from Equation 9.35 (or couple efficiency from Equation 9.42). In Table 9.1, the 
highest generator efficiency was found when « h = 3.62 V . 

Table 9.1 and Table 9.2 demonstrate the process both with and without interface resistances. Table 9.1 
demonstrates the calculation of a segmented p-type element without contacts, while Table 9.2 
demonstrates an n-type element with metal interconnects and interface resistances. For simplicity, the 
metal contact and contact resistances were given a 1 K budget (each) for the temperature drop. Such 
contacts have little effect on the value for u within the thermoelectric materials. The value for u in the 
metals, however, is temporarily reduced due to the Peltier effect at the interfaces (Aw = u 2 TAa). 

In order to calculate I and further operating conditions, the total heat flux UtotiilV-^total or P ower / area 
desired W/A tota j must be given. 

The current density J v is calculated from Equation 9.46, Equation 9.44, and Equation 9.49, giving 

1+ -2. 

T _ t/total,h A V , Q w > 

7 P ~ ~A X * ~®~ < 9 - 56 ) 

A total *p,h *n,h 

The n-element current density, / n , and I can then be calculated from Equation 9.45 and Equation 9.37. 

For example, assume we desire to make a thermoelectric couple out of the p-type element of Table 9.1 

with the n-element of Table 9.2. The area ratio which ensures that the desired u current flows in each 

\ 61.63 

element using the same electric current /, is from Equation 9.45: — = = 1.57. If the desired total heat 

& H A n 39.34 

1 
\y i + TT^ a 

flux into the hot-side is 20 W/cm , then from Equation 9.56, J„ = 20 T X ^^ = 34 =-. 

H p cm 2 0.4283V + 0. 5388V cm 2 
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Similarly, J n — — 53 — j. The desired length (from Equation 9.37) is 1= 



Th uKdT 39 - 34 — 
t c UKal _ cm . 



34- 



A 



1.16cm. 



For a total thermoelectric couple area of 1 cm ,A p +A n = 1.0cm 2 , the individual areas are A n = 0.39 cm 2 , 
A p — 0.61 cm 2 , and finally the current in the couple is I = 20 A (from I = JA). The output voltage for each 
couple is 0.142 V using Equation 9.43. 



9.5 Thermoelectric Compatibility 

If the compatibility factor s (most efficient u current, Equation 9.18) of one part of the thermoelectric is 

significantly different from the s of another part, there will be no suitable current where both parts are 

operating close to maximum efficiency. This is the physical basis for thermoelectric compatibility, and is 

most apparent for segmented generators. 

To achieve high efficiency, segmented generators use large temperature differences to increase the 

AT 
Carnot efficiency r/ c = . Since the material thermoelectric properties (a, k, p) vary with 

Th 
temperature, it is not desirable or even possible (most materials have a maximum operating temperature 
where they may melt or otherwise decompose) to use the same material throughout an entire, large 
temperature drop. Ideally, different materials can be combined such that a material with high efficiency at 
high temperature is segmented (Figure 9.6) with a different material with high efficiency at low 
temperature. 2 In this way both materials are operating only in their most efficient temperature range. 

If u could be constrained to be always equal to s, then the most efficient material to choose for a 
segment would be that with the highest thermoelectric figure-of-merit z. In this case, known as infinite 
staging (or upper limit of efficiency ), the interface temperature between segments would ideally be the 
temperature where the z of both materials cross. For example, according to Figure 9.8, the best infinitely 
staged p-leg (0°C to 1000°C) would contain (Bi,Sb) 2 Te 3 , Zn 4 Sb 3 , TAGS, CeFe 4 Sb 12 , and SiGe with 
interfaces of about 200°C, 400°C, 550°C, and 700°C. 

Unfortunately, in a real generator, u = s is not possible, so a compromise value for u must be selected. 
If the compatibility factors s of the segmented materials differ substantially, all segments cannot be 
simultaneously operating efficiently, and the overall efficiency may actually decrease as compared to a 
single segment alone. Figure 9.7 shows graphically that a suitable average value for u can be found for the 
three materials (Bi,Sb) 2 Te 3 , Zn 4 Sb 3 , and CeFe 4 Sb 12 , which have compatibility factors within about a 



Segmented Generator 



Cascaded Generator 



t5vWdR L 



^H 




FIGURE 9.6 Schematic diagram comparing segmented and cascaded thermoelectric generators. The cascaded 
generator has a cascading ratio of three. 
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FIGURE 9.7 Comparison of reduced efficiency as relative current density, u, varies for different p-type 
thermoelectric materials. An average value for u can be found for (Bi,Sb) 2 Te 3 ( 125°C), Zn 4 Sb 3 (300°C), and CeFe 4 Sb 12 
(550°C), which gives a reduced efficiency (indicated with a "•") near the maximum efficiency. SiGe (800°C), on the 
other hand, has such a low compatibility factor s that using a u appropriate for the other materials would result in a 
negative reduced efficiency for SiGe. This makes SiGe incompatible for segmentation with the other thermoelectric 
materials. 

factor of two. The reduced efficiency at this average u is not far from the maximum reduced efficiency. 
SiGe, on the other hand, has a much lower value for s, such that if the u shown in Figure 9.7 is used, a 
large negative efficiency will result for the SiGe segment and the overall efficiency will decrease. If a 
smaller u is used, so that positive efficiency will result from the SiGe segment, the efficiency of the other 
segments will have deteriorated more than the efficiency increase from the SiGe segment. Thus, despite 
having a reasonably high value of z for good efficiency, SiGe cannot be segmented with the other 
materials in Figure 9.7 because of different compatibility factors. 

As a rule of thumb, the compatibility factors of segmented materials should be within about a factor of 
two. Within this range, a suitable average u can be used, which will allow an efficiency close to that 
determined by z. Outside this range of s are materials which are incompatible where the efficiency will be 
substantially less than that expected from z. The compatibility factor is therefore, like z, a thermoelectric 
property essential for designing an efficient segmented thermoelectric device. 



9.5.1 Materials Selection 

For segmented generators high z materials need to be selected, which have similar compatibility factors, s. 
Other factors (not considered here) may also affect the selection such as: thermal and chemical stability, 
heat losses, coefficient of thermal expansion, processing requirements, availability, and cost. 

The compatibility factor (Figures 9.8 and 9.9) can be used to explain why segmentation of 
(AgSbTe 2 )o.i5 (GeTe) .85 (TAGS) with SnTe or PbTe has produced little extra power, but using filled 
skutterudite would increase the efficiency from 10.5% to 13.6%. 

Very high efficiency segmented generators to 1000°C could be designed with skutterudites or 
PbTe/TAGS as long as compatible, high-temperature materials are used. The compatible, high zT 
n-type material La 2 Te 3 4 would be ideal as long as a compatible p-type material is found. 

For the high-temperature p-type element, a high zT material that is also compatible with PbTe, TAGS, 
or Skutterudite, has not been identified. Even if the material has low zT, e.g., zT = 0.5, it will produce 
some power, as long as it is compatible. For a material with a low zT to be compatible with PbTe, TAGS, or 
skutterudite, it must have s > 1.5 V , ideally s ~ 3 V . Since s ~ z/2a, the zT ~ 0.5 material cannot 
be a high Seebeck coefficient band or polaron semiconductor. Materials with high z and s have 
thermoelectric properties typical of high a metals. In a metal, the thermal conductivity is dominated by 
the electronic contribution given by the Wiedemann -Franz law K e = LT/p where L ~ 2.4 X 10~ 8 V~/K . 
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FIGURE 9.8 (a) Figure-of-merit (zT) and (b) compatibility factor (5) for p-type materials. 

The compatibility factor s ~ a/(2/cp) ~ a/(2LT) would then be appropriate if a is greater than 100 /jM/K 
at 1000 K. For example, a candidate for such a refractory p-type metal is Cu 4 Mo 6 Se g . 

9.5.2 Cascaded Generators 

Cascaded generators (Figure 9.6) can avoid the compatibility problem between segments. In a 
segmented element all segments are thermally and electrically in series so that a single M h defines u( T) 
throughout the element. However, a cascaded device contains an independent electric circuit for each 
stage, allowing an independent / and therefore u in each stage. In this way, different optimal values of u 
can be used for each stage. The compatibility difference between materials makes cascading always 
more efficient than segmenting. ' However, cascading is much more difficult to implement than 
segmentation. 

Truly independent circuits will require an electric connector between a stage at high temperature and 
the load at ambient temperature. In practice, it is best not to connect the high-temperature stages directly 
to the load. Such connectors must have some loss because they can neither have very low electric 
resistance (or they will conduct heat away from the hot-side due to Wiedemann Franz law), nor high 
electric resistance (large Joule losses). It can be shown that the loss from such connectors is proportional 
to 1/N where N is the number of couples for each connector. Thus, to minimize these losses, the ratio of 
the number of thermoelectric couples to the number of connectors should be large. 

To avoid such losses entirely the electric current should pass from the high-temperature stage to 
the load by going through the thermoelectric elements of the low-temperature stage (Figure 9.6). 
The differing values of u is provided by having a different number of couples in each stage. The ratio of 
the number of couples (Nj/N^ in the cooler stage (Stage 2, N 2 couples) to that of the hotter stage (Stage 
1, Nj couples), called the cascading ratio, 2 is derived by equating the heat flux out of the hot stage to the 
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FIGURE 9.9 (a) Figure-of-merit (zT) and (b) compatibility factor (s) for n-type materials. 

heat input to the cold stage (thus T in the following formulas is the interface temperature between the 
two stages). This can be expressed succinctly (including both the Peltier and conduction terms) using the 
thermoelectric potential 23 Equation 9.49. 

Equating the heat from Ni stage 1 couples with N 2 stage 2 couples operating with the same electric 
current /, gives 






^P.l " n,l 



a l,p r - 



1 
s l,p 



<*l,nT 



1 

s l,n 



<b , 



± n 



1 1 

«2. P T H a 2 ^T - 

s 2,p s 2,n 



(9.57) 



where, for maximum efficiency, u is approximately equal to an averaged s. 

The simplification where the thermoelectric properties (a, p, k) are constant with respect to 
temperature is given by Harman. For SiGe cascaded with TAGS/PbTe the optimum cascading ratio is 
2.21, so that there should be about twice as many TAGS/PbTe couples as SiGe couples. With this 
cascading ratio the efficiency of the generator can achieve the sum of the two stages, 
10.39% + 4.39% = 14.78%. 24 

9.5.3 Functionally Graded 

The importance of compatibility has been made apparent for a segmented thermoelectric generator, but 
compatibility is also a consideration for all thermoelectric devices whether the materials are 
homogeneous or inhomogeneous from segmentation or being functionally graded. Since the 
compatibility factor of even a homogeneous material is temperature dependent, the change in s from 
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one end of the leg to the other will, in general, adversely effect device performance. The consideration of 
compatibility within the same material has been called self-compatibility. Because of the explicit 
temperature dependence of s (Equation 9.18) in conjunction with the typical temperature dependence of 
z and a, the problems with self-compatibility are generally more apparent at low temperatures. 

For inhomogeneous materials, the compatibility factor must be considered when optimizing the 
compositional variation. The thermoelectric properties of most materials can be altered with small 
changes in doping, which allows a tuning of both z and s. The most efficient material will not only have 
high zT, but will also be compatible with the other materials. For example, (Bi,Sb) 2 Te 3 with the highest 
zT also has a large compatibility factor compared to TAGS. By lowering the carrier concentration from 
the highest zT composition, a more compatible material can be used which is more efficient when 
segmented with TAGS despite having a slightly lower zT. 

True functionally graded materials can have a continuum of different compositions along the length of 
the material, which can enhance efficiency. ' Ideally, it would be best to have a local criterion for 
maximizing efficiency. By examining the reduced efficiency (Equation 9.15) at every location, it can be 
rationally determined if alterations in materials properties which change z and s will further increase 
efficiency. This also avoids the complex numerical methods used to predict the performance gains of 
functionally graded designs. 

9.6 Design Optimization 

The performance of a thermoelectric generator is dependent on many variables that could be optimized 
globally to find the optimum design. However, by using a reduced variable approach to the design 
problem, many interdependencies of the design variables are eliminated, which allows a better 
understanding of the effect of each variable. 

The following discussion is intended primarily for ideal thermoelectric generators. In "ideal" systems 
there are no thermal losses, no contact (thermal or electric) resistances, and only one-dimensional heat 
flow. Many nonidealities can be considered part of the system design, allowing the thermoelectric portion 
of the generator to be considered ideal. In addition, the effect of some nonidealities scale with length the 
same way as the thermoelectric material (such as interconnect metals with a fixed AT budget), which 
allows them to be incorporated without significant modification of the design algorithm. Finally, 
nonidealities which scale with length differently (such as contact resistances) should lead only to small 
corrections that will give a quick convergence if calculated iteratively. 

9.6.1 Efficiency Matrix 

The first goal of the design process is to evaluate the highest possible thermoelectric efficiency for all hot- 
and cold-side temperatures (of the thermoelectric generator, not the heatsinks), which may be viable. 
This will produce an optimized efficiency that is only a function of the thermoelectric hot- and cold-side 
temperatures: 

V = WOh, T c ) (9.58) 

The presumption is that any other variables (such as materials chosen, interface temperatures, geometry, 
current, etc.) that may be required for the calculation of efficiency can be optimized given a T h and T c . 
This is true for the thermoelectric material interface temperatures, but less true for size of metal 
interconnect and contact resistance. 

9.6.1.1 Interface Temperatures 

Once the materials are selected based on z and s, the interface temperature between thermoelectric 
segments must be optimized. The interface temperature can be estimated by where the value of z crosses 
or by using the approximation that u remains constant within a thermoelectric leg. 1 Calculating the 
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efficiency (as described above in Section 9.4.2) for different initial conditions M h for both the n- and 
p-element will allow a quick convergence to the most efficient Equation 9.42 u p >h and u n ^. 

The ideal optimum interface temperature is the one which has equal reduced efficiency for either 
material at that temperature (see 209°C rows in Table 9.1). By estimating the change in reduced efficiency 
with temperature (using an optimized u), the optimum interface temperature can be found with an 
iterative process. Often the interface temperature is simply the maximum operating temperature of a 
material (see 400°C rows in Table 9.1). 



9.6.2 Thermoelectric Element Length 

Once the optimized efficiency (Equation 9.58) is found given a T^ and T a the values of m(T) and <P(T) for 
both the n- and p-element are defined. Most of the remaining performance parameters also require the 
thermoelectric element length. This is usually determined by the desired total heat flux U totl i^/A tobl i (or 
power/area W/A tot3l ). As shown in the discussion of Equation 9.56, I can be calculated from UtotalV-^tatai 
with u(T) defined by optimizing the efficiency between T^ and T c . Thus, / is a function only of T^, T c , and 

^total,h'Aotal : 

I = l(T h , T c , C/ total /A total ) (9.59) 

Once functions Equation 9.58 and Equation 9.59 are evaluated for a variety of T h , T c , and l/ tota ij,/A tota i, 
they can be incorporated into the system model to find the optimal system operation condition. 

For example, often in a thermoelectric generator the power/mass is the primary concern. In this case, 
the power/mass can be increased by reducing the mass of the heat exchangers at the cost of reducing the 
temperature difference, which lowers the efficiency. By knowing how the optimum efficiency and length 
vary with input temperature and heat flux, the exact system solution is found without requiring the 
systems analysis to be capable of thermoelectric calculations. The system power and voltage are directly 
proportional to the size of the generator (through Aotai) and number of couples. Once the system trades 
are complete, the final configuration of the thermoelectric generator can be established. 

9.6.3 Voltage 

The voltage produced Vsystem is the number of couples connected in series N seT i es times the couple voltage 
Vcoupie (Equation 9.43): 

^system ^couple-^ series \J.O\J) 

Thus, the number of couples in series is determined by the voltage requirement. Often redundancy is 
desired by including additional parallel circuits N para n e i: 

-^system -^series-^parallel y^-®*) 

Once the thermoelectric length is fixed, the total power desired Wwill define the total cross-sectional area 
Aotal- The relationship between the area of a couple and the number of couples N system , is given by 

W 

^couple — TT (9-62) 

v L'total 



a '1^ system 

^total 

9.6.4 Maximum Power Density/Matched Load 

For many power generation applications, the power output is of greater concern than efficiency. In solid- 
state thermal-to-electric conversion, an important characteristic of the system is the aerial power density 
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(power provided per cross-sectional area W/A tota |). Equation 9.47 and Equation 9.50 show that the aerial 

power density is inversely proportional to the thermoelement length /. At a local level, this can be seen 

AT W 

using the approximations VT ~ — — and — ~ Pi. Equation 9.10 then becomes 

/ J\ 

W AT 2 

— ~ — : — X Ku(a — upK) (9.63) 

J\ I 

Thus, any aerial power density can be achieved by adjusting I, allowing the reduced current u which 
provides the maximum efficiency to always be used. 

In real generation systems it is the characteristics of the heat source and sink, heat exchangers, and heat 
concentrators that are needed to find the conditions for maximum aerial power density. Such passive 
heat transport systems require significant temperature drops to achieve high heat fluxes needed for high 
aerial power densities. The temperature drop used in the thermal transport system reduces the 
temperature drop in the thermoelectric generator, and therefore the Carnot efficiency. Thus, aerial power 
density of the system can be increased at the cost of reduced system efficiency. As a rule of thumb, the 
maximum power is generally found when only 1/2 of the heat source-heatsink temperature difference is 
used across the thermoelectric (by thermal resistance load matching ) with the other half used in the 
heat transport system. Nevertheless, the thermoelectric generator itself should be designed to operate as 
efficiently as possible with the heat flux and temperature difference it is allotted. In the discussions in this 
chapter, (AT) refers to the temperature difference across the thermoelectric only, not the entire system 
(see Chapter 11). 

Maximum power density in thermoelectric generators has been frequently analyzed ' ' by 
distinguishing the different operation condition which in a given thermoelectric generator provides 
more power (at lower efficiency) than the maximum efficiency condition (u = s). For a given generator 

(Z is fixed) and constant temperature difference (AT), Equation 9.63 the maximum aerial power density is 

Z 
found when u = — which is slightly larger than u = s (Figure 9.2). However, the larger u will require 
2a 

additional heat (from the Peltier effect) to be supplied and the overall process is at a lower efficiency 

(because of Joule losses). 

This analysis is traditionally performed using load resistance i^ad as the current adjusting variable as 

opposed to u or even electric current I. In addition, the approximation a, k, p, and therefore z are 

AT 
constant is also commonly made. By approximating VT ~ — — , and combining with /(-RLoad + R) = aAT 

and Equation 9.9 gives 



z R 
aRhozd+R 



(9.64) 



z 
Thus, the u = — condition is given by -RLoad — R (matched load) while the u — s condition corresponds 
2a 

to -RLoad =RyT+ZT (with Z — z), the load resistance for maximum efficiency. 

However, when designing a system (with length, I, as a design parameter), the matched load condition 

(J? L oad = R, u= — ) is a poor choice not only when considering efficiency, but also when optimizing 
2a) 

aerial power density, size, weight, or even voltage. For example, consider a design where the matched load 

condition was selected for the thermoelectric generator. This generator has Q heat flux supplied to it, at 

W . I z \ 

some T h and T c , and generates — = r)j_Q aerial power density 1 17 j_ is the efficiency when u = — ), 

A 2a \ la 2a) 

with a length lj_. Now, replace this u = — generator with a u = s(-R Loa d — R-Jl+ZT) generator 
2a 2a 
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having the same T h , T c , and heat flux Q. For a small AT generator, Equation 9.38 can be approximated 

kAT 

Q^ ——(aTu+l) (9.65) 



To utilize the same heat flux Q, the u = s generator will need to have a smaller length l s , than a 

z 

2a 
given by 



z z 

t( = — generator because s < — . From Equation 9.65 the ratio of the two lengths is approximately 
2a 2a 



L VzT+1 

if ~ ^T— ^ 

— + l 

Operating at the same temperatures but different currents, the longer, matched load generator will always 
have lower efficiency than that of the shorter u = s generator, because t/_z_ < t/ 5 . For comparison, a 

2a 

z . . zT 

generator operating with u = — will have a reduced efficiency of T/ r — , which is always less than 

2a 2zT + 4 

that of the u = s generator given by Equation 9.20. Since the heat entering both generators is the same, the 

( w ■ \ 

u — s generator will provide more power I — = 17Q, Q is the same ] as well as greater efficiency (and 

higher voltage). Also because l s < l_z_ the u = s generator will be smaller and therefore lighter. 

2a 

It is worth emphasizing again that, for a system which includes heat exchangers, the temperatures 
chosen for the thermoelectric will be different for the maximum power density and the maximum 
efficiency problem. However, in both cases, the generator should be designed such that the current used 
corresponds to the maximum efficiency and not the matched load condition. 

When losses are introduced in the thermoelectric portion of the generator which do not scale with 
length in the same way as the thermoelectric materials, the optimal current may deviate toward the 
matched load condition. An example is electric contact resistance which is independent of length. For 
systems with small /, the contact resistance becomes more important and so designing a longer generator 
by increasing the current toward the matched load condition increases both thermoelectric generator 
efficiency as well as power density. Such losses also lower the effective Z for the generator (an average of 
the z(T) of each material ) including contact resistances, which also brings the load for maximum 
efficiency closer to the matched load condition. 



9.6.5 Nonoptimal Operating Conditions 

Once the optimal configuration is established, the performance at nonoptimal conditions, such as the full 
I—V curve, can be calculated. In general, the heat flow, the temperatures, and the current can all vary 
from the optimal but the geometry remains fixed. The heat flow and temperatures will change in a 
correlated way determined by the thermal impedance of the components external to the thermoelectric 
converter, just as the electric current will change due to a change in the external electric impedance. Given 
the electric current in the generator and two out of three of the heat flux, hot- and cold-side temperatures 
(or equivalent relationships), the relative current density u, and therefore the generator characteristics 
can be determined. 

For example, if the hot- and cold-side temperatures are known (e.g., remain constant for low external 
thermal impedance) the relative current density u(T) of each element can be calculated (and estimated 
with ~) from 



— = UK I 

A )t c 



1-7=1 MKdT« ukAT (9.67) 
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The approximation in Equation 9.67 is good when uCT) remains approximately constant. This is typically 
true for electric currents up to those required for maximum efficiency. For very large currents, u no 
longer remains constant and diverges at the hot-side. A w h = ±00 is the current for a thermoelectric 
cooler operating at maximum AT (VT = at Th). The corresponding u c does not diverge, so it is more 
stable to compute using u h rather than u c as initial conditions. 

If the effective external thermal impedance is high, an I— V curve can be calculated assuming the heat 
supplied remains constant (and the hot- or cold-side temperature or a relationship is given). In this case, 
the three unknowns (the unknown temperature, and an initial value for u n and M p ) are solved from three 
equations: two of the form of Equation 9.67, for n-type and p-type, and the heat-flow Equation 9.49 for 
the couple. 

For low current operation, the Peltier cooling at the hot end of the thermoelectric will decrease, 
requiring an increase in the hot-side temperature and the interface temperatures. Such high-temperature 
operation may advance the degradation of the thermoelectric materials. If such degradation is 
detrimental to system performance, lower optimal operation temperatures should be selected. 

Because the geometry of the elements (in particular the length and area of each segment) remains the 
same when changing the current, the interface temperatures between the segments will drift from their 
optimal values. The interface temperatures can be found by finding the interface temperatures that keep 
the lengths of each segment constant, using Equation 9.67 as a guide. Even for low external thermal 
impedance, where the hot- and cold-side temperatures remain constant, the interface temperatures 
between the segments will change somewhat with varying electric current. 
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